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ABSTRACT 

We present a measurement of the lifetime of ground state atomic carbon, C( 3 P), 
against ionization processes in interplanetary space and compare it to the lifetime ex- 
pected from the dominant physical processes likely to occur in this medium. Our mea- 
surement is based on analysis of a far ultraviolet (FUV) image of comet C/2004 Q2 
(Machholz) recorded by the Galaxy Evolution Explorer (GALEX) on 2005 March 1. 
The bright C I 1561 A and 1657 A multiplets dominate the GALEX FUV band. We 
used the image to create high signal-to-noise ratio radial profiles that extended beyond 
lxl0 6 km from the comet nucleus. Our measurements yielded a total carbon lifetime 
of 7.1 - 9.6xl0 5 s (ionization rate of 1.0 - 1.4xl0 _6 s _1 ) when scaled to 1 AU. This 
compares favorably to calculations assuming solar photoionization, solar wind proton 
change exchange and solar wind electron impact ionization are the dominant processes 
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occurring in this medium and that comet Machholz was embedded in the slow solar 
wind. The shape of the C I profiles inside 3xl0 5 km suggests that either the CO life- 
time is shorter than previously thought and/or a shorter- lived carbon-bearing parent 
molecule, such as CH4 is providing the majority of the carbon in this region of the coma 
of comet Machholz. 

Subject headings: atomic data - atomic processes - comets: individual (C/2004 Q2 
(Machholz)) - interplanetary medium - methods: data analysis - molecular data - 
molecular processes -planets and satellites: auroras - Sun: solar wind 



1. Introduction 



The volatile content of comets makes their compositional study via remote sensing particularly 
interesting and effective. Dozens of molecules a nd their elemental constituent s have been detected i n 
cometary coma and ion tail emissions (see e.g.. IBockelee-Morvan et al.ll2004l ; iFeldman et al.ll2004l ). 
Determining cometary nuclear composition from these remote sensing studies requires detailed 
understanding of the physical processes that govern coma and ion tail emissions. The study we 
report on here concerns one of these physical processes: the ionization of ground state atomic 
carbon, C( 3 P). 

UV emission from atomic car bon in comets was fir s t detected in sou nding rocket observations 



of comet C/1973 El (Kohoutek; IFeldman et al 



1974 



Opal et al.l Il974l ) . The emission consists 



of resonant fluorescence fro m C( 3 P), with m ultiplets at 1561 A and 1657 A (see e.g., Figure 2 of 



Festoul Il984l or Fi gure 1 of 



Tozzi et al 



1998). These emissions are extremely extended spatially 



(e.g., Figure 1 of iMcPhate et al.lll999l ). implying a very long lifetime of C I against ionization 
processes. The distribution of the C I emission rules out a direct nuclear s ource. For most comets. 



CO a nd CO2 carry the bulk of carbon into cometary comas (e.g., Table 1 of IBockelee-Morvan et al 



2004J). CH 4 , C 2 H 2 , C 2 H 6 , CH3OH, and H 2 CO also contribute. As we discuss in detail in gj the 
final radial velocity of a carbon atom liberated from any of these molecules is ~4kms _1 . The 
lifetimes of the parent molecules are long enough to place the liberated carbon well beyond the 
collision sphere of any comet observed in the modern era. Thus, the average carbon atom released 
from a typical comet becomes an excellent test particle to probe conditions in interplanetary space 
unperturbed by the presence of the comet or any other solar system body. By observing the radial 
profiles of C I in comets and fitting them with coma models, we measure the lifetime of atomic 
carbon against ionization processes in interplanetary space. 

The majority of the C I measurements to date have been conducted by instruments on sounding 
rockets, IUE and HST, which have small fields of view (FOVs) relative to the typical angular extent 
of the C I emission. Exc eptions to this have been objective grating spectra of comet C/1975 Vl-A 
(West ISmith et al.lll98(lh. narrowband wide-fie ld images of comet C/1995 01 (Hale-Bopp) recorded 
by Polar spacecraft (jBrittnacher et al.ll200ll ). Wide Field Imaging Survey Polarimeter (WISP) 
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sounding rocket observations of comet Hale-Bopp ([Harris et al, 1 119991 ) . and imaging/spectro- imaging 
observations of Hale-Bop p by MSX (jKuppermanl 11999) . With a technique similar to the one 
we use here, iKuppermanl found that the lifetime of atomic carb o n was consistent with the solar 
photoionization lifetime of carbon calculated by iHuebner et al.l (|1992l ). However, comet Hale- 
Bopp was a comet of unusual size . Several physical proce sses, including collisional heating in 
the dense inner coma ( Combi et al. 1999 : [Harris et al. 2002 ) and gas phase chemical interactions 
between coma species (e.g., iGlinski et al.1 12004 ) acted to extend the distribution of coma species 
beyond what was expected using coma models based on comets with more typical production rates. 



Kuppermanl did not consider these effects. We use observations of comet C/2004 Q2 (Machholz), 



which is a much more "typical" comet. These observations suggest that the lifetime of carbo n 
against ionization is about half of the photoionization lifetime calculated bv lHuebner et al .1 ([1992). 
This implies that another process, which we show to be solar wind ionization, is comparable. 



2. Observations 

The Galaxy Evolution Explorer (GALEX) is a NASA Small Explorer (SMEX) mission designe d 



to map the history of star formation in the universe ([Martin et al.ll2005l : lMorrissey et al.ll2005l . 120071 ) . 
Its 1.2° FOV diameter and high sensitivity to extended sources also makes GALEX well suited to 
cometary coma studies. GALEX operates in two modes: direct imaging and objective grism 
imaging. Using a beam splitter and two detectors, GALEX images are simultaneously recorded in 
two bands: 1350-1750 A (FUV) and 1750-3100 A (NUV). 

GALEX made both direct-mode and grism-mode observations of comet C/2004 Q2 (Mach- 
holz). The comet parameters are shown in Table [TJ Over the course of the observations, the comet 
moved a substantial fraction of the GALEX 1.2° FOV, so we quote only the center of the spacecraft 
FOV, which was used for all observations. A log of observations is given in Table EJ which includes 
a background direct-mode image that was recorded a month earlier. Pipeline processing did not 
complete properly for grism orbit 3 due to a telemetry dropout, so we excluded it from our analyses. 

The top two panels of Figured] show images created by the GALEX pipeline processing system. 
The 1657 A multiplet is the dominant feature at all distances from the nucleus. Stars that appear 
as point sources in the direct-mode image are streaks in the grism-mode image since each stellar 
spectrum formed by the objective grism is individually imaged on the detector. On this scale, 
the comet motion is not readily apparent, but must be considered in order to provide high-quality 
radial profiles of the direct-mode image and co-add all of the grism images. 



3. Data Reduction 



The GALEX MAMA imagers detect and time tag individual photons. The photon lists are 
telemetered to the ground along with the spacecraft aspect information. Using reference stars in 
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Table 1. Observational Parameters 



Parameter 


Value 


Comet observed 


2005-03-01 


Heliocentric distance 


1.3 AU 


Heliocentric velocity 


11.0 kms" 1 


Heliographic lat. 


30.3° 


Heliographic Ion. 


141.0° 


Geocentric distance 


0.75 AU 


FOV R.A. 


4hr 13'50," 


FOV Dec. 


82°45'10" 



Table 2. Log of GALEX comet C/2004 Q2 (Machholz) observations used in this work 



Type 


Name a 


Exp 
time (s) 


Date 


Start 
Time 


Back, direct 


1-fd-x.fits 


271 


2005-01-30 


03:04:42 


Comet direct 


2-fd-x.fits 


1635 


2005-03-01 


19:29:19 


Comet grism 


1-fg-x.fits 


1618 


2005-03-01 


11:16:19 


Comet grism 


2-fg-x.fits 


1619 


2005-03-01 


12:54:57 


Comet grism 


4-fg-x.fits 


1619 


2005-03-01 


16:12:13 


Comet grism 


5-fg-x.fits 


1619 


2005-03-01 


17:50:51 



a fllcnamc prefixes arc GI1.104001.COMET_2004Q2.000 




Fig. 1. — Top left: GALEX default pipeline processed direct-mode comet Machholz FUV image. 
Top right: the same, but grism-mode orbit 1. Bottom left (direct-mode) and right (grism-mode): 
images reconstructed in the frame of the comet, with background sources removed, and diffuse 
background subtracted, as described in Sj3j Four grism-mode observations are co-added (note 
trailing of FOV). The grism-mode images have been oriented so that the dispersion direction is 
horizontal and longer wavelengths are found to the right. Note that three bright stars off the edge 
of the FOV have been dispersed by the grism into the FOV. 
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the FOV, the astrometric solution is refined to better than 0.6" (jMorrissev et al.l 120071 ). In the 
reference frame of the stars, comets are trailed during the ~ 1500 s GALEX observing segments 
("visits"). When several orbits are available, comets appear multiple times across co-added images. 
To eliminate comet trailing and multiple comet images, we used astrometrically corrected photon 
lists and an accurate comet ephemeris to cast photons into a comet frame of reference, which trails 
the stars. We obtain high-precision comet ephemerides from the JPL HORIZONS system using 
"@galex" as the observatory, which ensures that effects such as GALEX orbit parallax are properly 
considered. All coordinate transformations were done using the IDL ASTRO routines Q 

Comet motion across the GALEX FOV and a mask that blocks out background sources allows 
us to completely remove flux from identified background sources. We used direct-mode observations 
of the same star field recorded a month before the Machholz arrived to create a high quality 



background source mask. Sources were identified by Emmanuel Bertin's SExtractor packaged that 
is run as part of the normal GALEX pipeline. 

For the direct-mode observation, masks were created using the extreme corners of the sources 
tabulated in the SExtractor catalog columns XMIN, XMAX, YMIN, and YMAX. The mask was 
used to blank out the background sources in separate count and exposure time images that were 
accumulated as the photon lists were read. The mask was applied every time the predicted position 
of the comet moved 1 pixel relative to the background sources. The GALEX flat-field image was 
applied to the exposure image in a similar way. When the end of the photon list was reached, the 
count image was divided by the exposure image to produce a flat-fielded count rate image with 
the background sources effectively erased. The FOV of this image was extended in the direction of 
comet motion in the same way that the stars were trailed when we did not mask them. 

For the grism-mode observations we converted the individual source masks into rectangles 
elongated in the dispersion direction and proceeded with the reduction in the same way as with the 
direct-mode for each visit individually. Photometrically correct addition of the four grism visits 
was achieved by adding count and exposure-time images and then dividing the results. Because the 
grism can disperse light from stars that are out of the FOV, our method of working from sources 
cataloged in the direct-mode image did not always produce masks for stars. These stars could 
easily be masked by hand or by modeling the comet emission (described in subtracting it from 
the comet image cast in the reference frame of the stars by a procedure similar to which the stars 
are masked and inspecting the residuals for sharp features. For the purposes of this work, stars left 
in the image did not prove to be problematic. 

In order to assure that our method of processing the GALEX photon list was consistent with 
that of the default GALEX pipeline, we ran our software without the comet ephemeris and without 
masking the stars. We found minimal count rate differences between our images and pipeline images, 



1 http://idlastro. gsfc.nasa.gov/ 
2 http://sextractor. sourceforge.net/ 
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except within 100 1.5" pixels of the edge of the circular FOV. We attribute this to details in the 
flatfielding process. In order to avoid dete ctor burn in, the GA LEX spacecraft executed a 1.5' 
spiral motion during pointed observations (jMorrissev et al.ll2005l ). The default GALEX pipeline 
applied the detector flat in the detector coordinate system every few seconds while images were 
being reconstructed in an astrometric coordinate system. The effect of the spiral motion on the flat- 
fielding at the edge of the FOV was therefore correctly considered. The flat-field provided to guest 
observers by the GALEX project is an average of the detector flat over each orbit. Our method 
of applying this flat therefore did not take into consideration the spiral motion. The flat changes 
slowly as a function of position except at the edges, so the differences in flat-fielding methods does 
not effect the quality of the data except at the edges. As the GALEX images are 3840x3840 pixels, 
we simply masked out the 100 pixel wide annulus when we were applying the flat-field corrections 
in our procedure. 

The lower-left panel of Figure [1] shows the resulting direct-mode image with the background 
discussed in the next section subtracted. The removal of discrete sources has been highly successful. 
The lower-right panel of Figure [1] shows the corresponding grism-mode image, which is a co-add of 
4 visits. 



3.1. Background 

Since the comet emission filled the field-of-view, it was not possible to estimate the contribution 
of the background from the comet image itself. We used the background image recorded a month 
before the comet image for this purpose. Discrete sources were removed using the mask generated 
from the SExtractor catalog. In order to fill in the missing data where the sources were removed, 
we smoothed the image using a 100x100 pixel boxcar average. The resulting image was translated 
to the average R.A. and DEC of the comet during the direct-mode observation. 

Before subtracting the image from the data, we checked for a ny effects that might contribute to 



variation of the background level over the course of one month. ISuiatha et al.l (|2009l ). in studying 
multiple GALEX observations of the Sandage nebulosity show that in the FUV, airglow from the 
1356 A and possibly 1304 A lines is responsible for the primary variation in the total GALEX count 
rate as a function of time over short and long timescales. Figure [2] shows the count rate in the 
GALEX FUV detector for the comet (upper curve) and background (lower curve) after removal 
of detector effects such as hot spots. The sinusoidal variation is induced by the combination of a 
bright star near the edge of the FOV and the spiral motio n of the GALEX spac ecraft. Otherwise, 



the count rate curves resemble the parabolas described by ISujatha et al.l (|2009l ). We adapted the 
their method to estimate the difference between the airglow contribution in the comet image and 
the background image: -5.66 x 10~ 5 counts s -1 pixel -1 . The primary contribution to the difference 
in the airglow contribution came from the fact that the background image was recorded only on 
the steep part of the count rate curve, so its average airglow was high compared to the comet 
observation, which covers the entire parabola. 
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Fig. 2. — Count rate in the GALEX FUV detector for the comet (upper data points) and back- 
ground (lower data points) visits. The count rate has been corrected for detector effects such as 
hot spots. The sinusoidal variation in the signals is due to a bright star near the edge of the FOV 
and the spiral motion of the GALEX spacecraft. The model described in the text plus the sinusoid 
is overplotted in red in both cases. The model with B = 1 (nightglow only) is overplotted to 
demonstrate the relative contributions of the day and nightglow signals. The average astrophysical 
signal without airglow is plotted as the straight lines for the two cases described in the text: po = 
0.7 (solid) and po = 0.5 (dashed). 
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We note that the ISujatha et al.l ()2009l ) method is largely empirical and must be repeated for 
each observing geometry. The comet and the Sandage nebulosity were both high-dec sources, so 
we expect our use of their parabolae to be reasonably accurate. However, the sensitively of our 
lifetime measurement to the precise value used for the background motivated us to investigate the 
airglow effect more carefully. We demonstrate here that a simple atmospheric model, geometric 
considerations of earth shadow on the atmosphere and accurate consideration of observing geometry, 
explains the variation in GALEX count rate seen during each eclipse period. Much like photometric 
data taken at various zenith distances, our method can be used to independently remove the effects 
of airglow from any GALEX eclipse. 

We constructed a model atmosphere using an exponential rooted in the F-region (altitude 
300 km) with an initial density po (detected GALEX FUV counts per km at an altitude of 300 km) 
and a scale height, h. We assumed that the distribution was the same over every point on the 
Earth but we defined a unitless "day boost" parameter, B, which multiplied the density of emitters 
when they were in sunlight. We used the orbit of GALEX described by two line elements (TLEs) 



archived at |http://www.spacetrack.org to determine the position of the satellite as a function of 
time, a solar ephemeris provided by the SUN-POS.PRO module of the IDLASTRO routines, and 
a simple stepwise integrator to compute the airglow brightness for each line of sight as a function 
of time into the visit. The final component of the model was the astrophysical count rate, which 
we assumed to be constant over the visit. The red solid lines in Figure [2] show the resulting best-fit 
model, including the sinusoid discussed above. Our model gave good fits for a range of parameters 
corresponding to airglow differences of -5.57xl0 -5 counts s _1 pix -1 to -3.6xl0 -5 counts s _1 pix -1 . 
An interesting extension to this exercise would be to use more sophisticated ionospheric and plas- 
maspheric models to derive the airglow emitter profiles and increase the accuracy of the absolute 
airglow determination. 

Having derived a best-fit background image, we subtracted it from the comet image. There 
was a slight oversubtraction around the edges due to the background image smoothing and the 
elongation of the comet image. The effected points were removed with an annular mask 175 pixels 
wide. Residual emission around the hot spot in the southwest portion of the image and a bright 
source toward the west were also masked. The resulting image is shown in the lower- left panel of 
Figure [TJ The co-added grism-mode images formed by our reduction procedure are shown in the 
lower-right panel of Figure [U with diffuse background determined in subtracted. 



4. Analysis 



We used the combination of GALEX direct and grism-mode data to measure the carbon 
ionization lifetime. We began by assuming that C I emission dominates the passband at all come - 
toc entric distan c es an d that CO is the sole parent molecule. As discussed bv lHuebner et al.l (|1992l ) 
and lTozzi et al.1 ( 19981 ) . CO dissociates in significant amounts into both ground state carbon, C( 3 P) 
and the metastable state C( X D). C( X D) decays to C( 3 P) in 4077s (a rate of 2.45xl0" 4 s" 1 ), inde- 
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pendent of heliocentric distance ( Hibbert et al.lll993). C( 1 D) can be phq toionized with a predicted 
rate that is ~10 times faster than that of C( 3 P) ( Huebner et al.lll992l ). However, the resulting 
photoionization lifetime, 2.8xl0 5 s (rate of 3.6x 10 -6 s _1 ), is still long (rate is slow) compared to 
the metastable lifetime (rate), so for the purposes of this work, we neglect the time spent in the 
metastable state and assume all of the CO dissociates directly to C( 3 P). In we explore the 
differences in the photochemical ejection velocities of C( 3 P) and C( 1 D), which could, in principle, 
affect our results, and find them to be small (~25%). 



We fit two-component lHaserl models to radial profiles formed from our background-subtracted 
direct-mode image (|4.14 Fig. |3|). We did this as a function of airglow offset in order to estimate the 
effect of our background subtraction uncertainty on the carbon lifetime (Fig. We find that the 
lifetime measured by this method agrees well with the theoretical ca rbon ionizati o n lifet ime calcu- 
lated using values for solar wind and photoionization tabulated by iRubin et al.l (J2009J) assuming 
comet Machholz was in the slow solar wind (Fig. HJ Table |3|) . 

Figure [3] shows radial profiles of four quadrants of the background-subtracted direct-mode 
comet image shown in the lower-left panel of Fig. [TJ Quadrant 4 is centered on the dust tail, 
quadrant 1 is oriented 90° counter-clockwise from quadrant 4, etc. T he profiles we r e crea ted using 
an adaptive binning ring-sum algorithm similar to that described by lHarris et al.l (|2002l ). For the 
profiles shown in Figure O at least 1000 counts per radial bin were used. For our fitting exercises, 
there were at least 10,000 counts per bin. The profiles show good agreement except for the one 
oriented toward the dust tail. 



4.1. The Haser model 



Haserl (jl957l ) first solved the simple differential equation which describes the radial distribution 



of a species traveling radially away from a comet nucleus at velocity v^ik and with lifetime against 
destruction, r„. A similar di fferent ial equation can be solved for the case of a daughter species, 



resulting in the 2-component Haser model (e.g.. iKrishna-Swamv &: Brandtl ll986). A critical com- 
ponent of the 2-compone nt lHaserl model is the ejection velocity, v e , of the da ughter from the parent 
molecule. iFestoul (|198ll ) pointed out that the common usage of v e in the lHaserl model was un- 
physical: the daughters were frequently assumed to be ejected radially outward from their parents, 
when in fact, they would be emi tted isotropically from the point of destruction of the parent. The 
vectorial model oflFestoul (119811 ) addresses this shortcoming. Many workers have noted that the 



2-component lHaserl model provides go od fits to comet profiles of daughter species. The problem 



is mapping the input parameters of the lHaserl model (e.g., parent and daughter outflow velocities) 



ity). 


Combi et al. 


Hasei model to be 



formula? in all our calculations. 
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Fig. 3. — Quadrant-by-quadrant radial profiles of comet Machholz FUV emission. The left-hand 
Y-axis shows the background-subtracted average count rate per pixel in each radial bin. The right- 
hand Y-axis shows the count rate converted into C I column density using the g-factors described 
in the text and the assumption that all of the emission in the FUV band is from C I. Quadrant 4 is 
centered on the antisunward or "Dust tail" direction, as indicated in the lower-left panel of Fig. [TJ 
Quadrant 1 is counterclockwise from quadrant 4, etc. Statistical error bars (total vertical extent 
2-0") are shown and are comparable to the line widths. 
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As outlined above, our starting point was to assume CO is the sole parent of C I. We also as- 
sume CO is produced within the collision sphere and thermalized to the bulk outflow velocity, t>buik> 
of the coma. Preliminary analysis of our contemporaneously recorded GALEX NUV data suggests 
that comet Machholz had a water production r ate i n the log(Q(HoO)) = 28 - 29 range. According 
to the formulae given by iBudzien et al.l (|1994l ) and iTseng et al.l (|2007l ) , this translates into t>buik 
= 0.74 km s -1 and 1 .02 km s -1 , respectively at the c omet's heliocentric position. Combining the 
Huebner et al. ( 1992 ) CO photodissociation rate and Rubin et al. (j2009l ) CO solar win d ionization 



rates results in a total solar quiet CO lifetime of 7.9xl0 5 s (rate of 1.3x 10~ 6 s" 1 ). Ilbzzi et al 



(j 19981 ) use an improved set of calculations and find a solar maximum CO photodissociation life- 



time of 4xl0 5 s (2.5x10 6 s 1 ). We estimate the solar minimum CO lifetime to be < 5xl0 5 s 
(> 2xl0- 6 s- 1 ). 

As discussed above, we can neglect ionization out of the C( X D) state because the metastable 
lifetime is so much shorter than the ionization lifetime, however, we must check to see if the 
ejection veloci ties, of C(3P ) and C( X D) from CO are sufficiently different to re quire indivi dual 
consideration. iHuebner et al.l calculate v P . = 4.9kms _1 and 3.9 k m s -1 , resp ectively. iTozzi et alJ find 



the ejection velocity, v e , of C(3P) from CO to be 4kms _1 . ITozzi et al.l do not directly calculate 
v e values for C( X D), but rather adopt a range of values from ftmik to 3.9 km s -1 . We used a less 
conservative lower l imit and assum ed that the C( X D) v e value is within 25% of the C(3P) v e value, 
as suggested by the IHuebner et al.l calculations. As discussed below and shown in Figure HI this is 
comparable to the uncertainty in our fitting and background removal activities, so we neglect it. 



The formula of IWhipple &: Huebnerl (|1976l ) shows that the collision sphere, over which gas 
interactions transfer excess energy between gas species, was ~2500km for Machholz. This is much 
smaller than scale lengths of CO or any of the major carbon parent molecules. Thus, the mo 
tion of carbon is d ominated by its ejection velocity. Since v e is significantly larger than i>b u ik) the 
Combi et al.l (|2004l ) formulas show that the daughter lifetime that would be der ived from sophisti- 
cated models is essenti ally unaltered from that which is derived from the lHaserl model. Conversion 
of our measured lHaserl scale lengths to the C I lifetime is a simple matter of dividing by the assumed 
ejection velocity. 



We used MPFIT ( More 



1978 



Markwardtl [2009) to fit two-component lHaserl models to the 



profiles. Because the statistical errors on the profiles we were fitting were so small, the formal 
errors of all the parameters were always less than 1%. The of f-diago nal elements of the correlation 
matrix were also always less than 1%. The two-component lHaserl model was able to accurately 
reproduce all of the profiles in Figure [3] except for the inner-most point, which consistently fell 
above the model. 

We experimented with excluding data below small cometocentric distances, r m ; n to see if this 
effected the C I lifetime. As r m j n was increased from to 1 x 10 5 km, there was a 20% decrease in the 
measured C I lifetime. Finding this to be an acceptable range of uncertainty in the fitting exercise, 
we choose the midpoint, 5xl0 4 km as our best value and proceeded with the other experiments 
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described below. 

Next, we varied r max , the upper cutoff of the data we fit, and the assumed airglow offset value. 
We created contour plots of the C I lifetime in the plane of airglow offset versus r max for a variety 
of combinations of quadrants. We found that quadrant 1 resulted in the most consistent behavior 
as a function of r max . For all 

r max. > 3x10^ km, in quadrant 1 all the measured C I lifetime values 
fell within 10% of each other. This is important. It means that the GALEX FOV, which extends 
to ~l.lxl0 6 km is large enough to reliably measure the carbon scale length, even though at this 
heliocentric distance, it was larger than the GALEX FOV. 

Figure |4] summarizes our fitting exercises. The C I lifetime is plotted against the assumed 
offset between the comet and background airglow levels. The plot was made by fitting quadrant 1 
data between r m j n = 5xl0 4 km and r max = l.lxl0 6 km (the largest radial bin containing at least 
10,000 counts). It is at the edge of the radial distribution that the background subtraction accuracy 
becomes particularly important, as any over /undersubtr action is directly reflected as a change in 
the number density of CI and therefore the scale length. We mark the two values of the airglow 
offset level discussed in £13.1l with the solid and dashed lines. The resulting 1 AU C I lifetimes are 7.1 
and 9.6xl0 5 s for the "best" and minimum cases, respectively. These correspond to total ionization 
rates of 1.4 and 1.0xl0 -6 s , respectively. Table [3] summarizes these results. 



4.2. Comparison to physical processes 



Photoionization is generally assumed to be the dominant ionization process for most cometary 
molecules and radicals, however, for long-lived species such as carbon and oxy gen, solar wind 
ionization processes become non-negligible (e.g., lAxfordl Il972l ; iRubin et al.ll2009i ). We find that 
for carbon, solar wind ionization processes are comparable to photoionization. Table [3] summarizes 
calculations of carbon ionization lifetimes (rates) from the dominant physical processes for a variety 
of solar conditions. We observed comet Machholz near the end of solar cycle 23, so we assume quiet 
sun conditions from the perspective of the photon field. The heliographic latitude of Machholz was 
30°, so Machholz was likely on th e upper edge of the slow solar wind zone ide ntified by Ulysses 
and SOHO SWAN observations (jSmith Sz Marsdenl Il995l : iBzowski et al.l 120031 ) . Our measured 



ionization lifetime of 7.1 - 9.6xl0 5 s (1.0 - 1.4xl0~ 6 s -1 ) therefore compares favorably with the 
predicted lifetime of 8.2xl0 5 s (1.2xl0~ 6 s -1 ). In Figured! we present a graphical interpretation of 
these results combined with the presentation of our most easily identified systematic uncertainty: 
determination of the relative contributions of airglow in the comet and background images. 



4.3. Grism Image analysis 



In £ j4.1l we described iHaserl model fits to the direct-mode profiles and showed that the measured 
C I lifetime was relatively stable as a function of r m i n and r max (20%). However, these fits all 
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Fig. 4. — Carbon lifetime at 1 AU vs. assumed background offset between the comet and background 
direct-mode images. The best-estimate value of our background offset calculated in M3.ll to be - 
5. 57x 10 -5 counts s -1 pixel -1 is shown as the solid vertical line. The resulting carbon lifetime is 
similarly indicated graphically and labeled "Best." The minimum background offset is indicated 
in this way with the dashed lines. The calculated values of the carbon lifetime from the first two 
columns of Table [3] are shown as dotted lines. The data indicate that during our observations, 
comet Machholz was embedded in the slow solar wind. 
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required a parent scale length that was extremely short, indicating the presence of contamination 
from the shorter-lived species in the GALEX FUV bandpass. We use a preliminary analysis of 
the GALEX grism-mode data to identify the species contributing the contaminating emission and 
demonstrate that this contamination is unlikely to effect our C I lifetime measurement. 

The general method we use for analyzing the grism-mode images is forward modeling. We 
create a radial profile of the expected emission in a particular line or spectroscopically unresolved 
band. Our software is str ucture d so that any model profile can be used, but for these analysis, we 



us e the 1 or 2-c o mpon ent iHaserl model, as appropriate, with scale lengths modified by the formulas 
of ICombi et al.l ([20041 ) to emulate the output of more sophisticated coma models. The profile is 
scaled by the product of the production rate, the g-factor of the feature in question (see below 
and Table HI) and the effective area of the GALEX instrument at the wavelength of the feature 



(|Morrissev et alJ[200 7^l. We then create a figure of rotation from the profile, thus transforming the 



profile into a coma image. We offset the image according to the GALEX grism dispersion relation 
from the astrometric center of the comet, determined by our reduction procedure described in Sj3j 
Next, we convolved the image by an estimate of the GALEX point spread function in grism-mode, 
which we create by making a figure of rotation from the cross-dispersion profile model generated 
from stellar spectra in the grism-mode image by the GALEX pipeline processing system. We then 
repeated the process of converting the model profile into a coma image for each of the three primary 
GALEX grism orders and summed the results. 

We calculated ^-factors at 1 AU and a heliocentric velocity of 11 kms -1 using a high-resolution 
solar spectrum from the Solar Maximum Mission normalized to absolution fluxes measured by 
SOLSTICE together with absorption oscillator strengths available from NIST. For the purposes 
of these calculations, we assumed optically thin conditions. The implications of this assumption, 
particularly for CO, are discussed below. The CO ^-factors were calculated assuming an excitation 
temperature of 70 K. The (^-factors and GALEX grism mode effective ctrccis (ire listed in Table HI 

We determined the background level in our grism-mode data by subtracting a preliminary 
model of the C I emission based on the direct-mode analysis. We assumed that CO was the sole 
parent of C I and used the same C I production rate, Q(C) = 2.8xl0 28 s~\ found in our direct-mode 
analysis. After experimenting with Q(C), and the C I scale length, we found that the grism-mode 
background was much higher than expected assuming the value from the direct-mode image and 
the ratio between the total efficiencies of the two modes. We attribute the high background to 
emission from the unmasked sources discussed in £j3j We found a best-fit background level of 
5.4x10 counts s _1 pixel -1 . The top two panels of Figure [5] shows the inner 2xl0 5 km of our 
GALEX grism-mode data before (left) and after (right) the C I subtraction. The positions of the 
major features in the primary GALEX FUV grism order have been indicated. 

After subtraction of the C I emission, a pseudo-continuum was evident (upper right-panel of 



3 http://galexgi. gsfc.nasa.gov/docs/galex/Documents/PostLaunchResponseCurveData. html 
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Table 3. Carbon Ionization Lifetimes (Rates) at 1 AU 







Ionization Lifetime x 10 5 s 


(Rate xl0~ 6 s -1 ) 




Process 




Quiet Sun 


Quiet Sun 


Active Sun 


Reference 






Slow Wind a 


Fast Wind b 


Slow Wind a 




C +hv -> 


C+ + e- 


24 (0.41) 


24 (0.41) 


10 (0.92) 


Huebner et al. ( 1992) 


C + H+ - 


> C++ H 


17 (0.59) 


40 (0.25) 


17 (0.59) 


Rubin et al. (2009) 


C +e" -> 


C+ + 2e" 


48 (0.21) 


20 (0.05) 


48 (0.21) 


Rubin et al. (2009) 


Total predicted 


8.2 (1.21) 


14 (0.71) 


5.8 (1.72) 




Measured 




7.1 - 9.6 (1.0 


-1.4) 







& v = 400 km s 1 , n e = 10 cm 3 
h v = 750 km s -1 , n e = 2.5 cm -3 



Table 4. g-i actors at 1 AU 



Feature 


g-factor a 


Eff. Area b 


C I 1561 A 


0.93 


18.9 


CI 1657 A 


4.22 


11.5 


CO total 


0.16 


14.9 


S I 1425 A 


0.29 


8.9 


S I 1475 A 


0.25 


21.6 


S I 1813 A 


4.90 


2.5 



a xl0 5 phot s 1 particle 1 
b Grism-mode peak order, cm 2 
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Fig. 5. — Machholz FUV grism-mode analysis. The frames show the inner 2 x 10 5 km of the co-added 
comet Machholz GALEX grism-mode images in various stages of residual subtraction. Background 
described in the text has been subtracted from all frames. The primary grism order positions of the 
brightest features in the FUV bandpass have been marked. Top left: no model subtraction. Top 
right: subtract C I from a CO parent with Q(C) = 2.8xl0 28 s~ 1 . Bottom left: also subtract CO, 
with the profile described in the text and Q(CO) = 3.0xl0 28 s _1 . Bottom right: modify original 
Q(C) to 1.8xl0 28 s-\ Q(CO) = 2.8xl0 28 s- 1 and subtract Q(C) = l^xlO^s" 1 from CH 4 . To 
increase the signal to noise ratio in the residual images, they have been binned by 3x3, 3x3, and 
4x4 pixels, respectively. Also note changes in color bar limits. 
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Fig. [5|). Since dust emission near the nucleus is negligible in the FUV, we took this to be emission 
from the more than 80 lines of the CO Fourth Positive system that fall in the GALEX FUV band- 



McPhate et al. 


1999: 


Lupu et al. 


2007 


)■ 


Lupu et al. 



Fourth Positive emission as a function of cometocentric distance to 3,000 km with and without 
considering self-absorption (optical depth) effects. They find evidence for detectable optical depth 
effects in the inner ~1000 km of comets with production rates similar to that of Machholz. Perhaps 



not surprisingly, when we subtracted a one-component iHaserl model of the CO emission from our 



C I residual, there was obvious oversubtraction at cometocentric distances <10,000km. We emu- 



lated optical depth effects by using a 2-component IHaserl model with an effective parent lifetime 



of 10,000 s (rate of lxlO^s" 1 ) at 1AU. The resulting residual, assuming Q(CO) = S.OxlO^s" 1 
is shown in the lower-left panel of Figure [5j Since we assumed optically thin conditions, our fit 
is not sensitive to whe ther or not CO comes from a nucleus- nat i ve source or is formed from a 
distrib uted source (e.g., |E berhar dt et al.lll987l ; iDi Santi et al.ll200ll : iBrooke et al.ll2003l ; iLupu et al 



2007: Bockelee-Morvan et al 



201CJ). We look forward to further work on this subject as we incorpo- 
rate radiative transfer calculations in our CO profile calculations and improve our grism background 
source cleaning technique. 

After subtracting the primary C I and CO sources, residual emission, centered roughly on C I 
1657 A was evident. This was not an artifact of the unexpectedly large grism background discussed 
above. It is reasonably well described by C I from a parent with an effective 1 AU lifetime of 
lxl0 5 s (rate of lxl0 _5 s _1 ), assuming CO has a lifetime of 5xl0 5 s (2xl0~ 6 s -1 ; see §4.1j) . We 
experimented with various CO lifetimes ranging from 2x 1 5 s (which happened to describe a large 
section of our data well) to 7.9xl0 5 s ( Rubin et al.ll2009l ). These CO photodissociation lifetimes 



correspond to rates of 5 and 1.3xl0 -6 s _1 , respectively. The problem with a very short CO lifetime 
is that if all of the C I came from CO, the model clearly over-subtracted inside of 2x10 km. From 
this we deduced that the dominant carbon parent molecule cannot have a lifetime as short as 
2xl0 5 s (5xl0" 6 s _1 ) and that at least one more carbon parent molecule was needed. 

Knowing that we needed at least one other carbon parent molecule besides CO to provide 
C I emission inside 2xl0 5 km, we t urned to infrared o bserv ations of comet Machho l z to p rovide 
suggestions. Based on the work of iBonev et al.1 (|2009l ) and iKobayashi &: Kawakital (|2009l ). CH4 
and CH3OH were the most likely candidates. Our estimates of the maximum production rates of 
atomi c carbon from CO and these molecules is shown in Table \E\ Of particu lar note is our adop tion 
of the iHuebner et al.l ( 19921) CO bran ching ratios even though we use the iTozzi et al. I { 19981 ) CO 
lifetime. According to Hueb ner et al.l . a significant amount of CO is lost to ionizatio n, leaving only 



50% for C I emission. Tozzi et al. 



find a significantly shorter total CO lifetime than IHuebner et al 



but do not calculate rates into the various branches. 

The d etailed photodi s sociat ive lifetimes and excess energies of the CH4 to carbon chain are 
not know. IHuebner et al.l dl992) calculate the photo rates of the various channels of CH4 to its 
daughters and the excess energies of these reactions. The total solar quiet lifetime of CH4 is 



Table 5. Carbon Parent Molecule Candidates 



Parent Q(C) max Q(C) r p v e 

Molecule (xlO 2 ^ 1 ) (xl0 5 s) (kms" 1 ) 



co a 


100 




CO 


50 b 


1.4 


CH 4 


35 d 


1.0 


CH3OH 


50 d 


1.4 



2.8 a 

1.8 <5 C 4 C 

1.2 e 1.4 >3.1 

>5.9 4.5 
2.8 f 8§ 



1 CO Fourth Positive emission 



TOuebner et al.1 $1992 ) 
c See UJ\ 



^Boney et al.1 (|2009h and iKobavashi Kawakital (|2009h 
e Decreases to 1.0xl0 28 s _1 assuming CO lifetime of 3xl0 5 s 
f From direct-mode profile fit 
g Carbon lifetime from this work 
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1.3xl0 5 s (7.7xl0 _6 s _1 ). The primary channel to at omic carbon is through C H2. The interstellar 



and solar photo rates of CH2 have been calculated by Ivan Dishoeck et al.l (|1996l ). who find the total 
lif etime is 7 - 11 x 1 3 s ( 9 - 14x 10~ 5 s -1 ). The majority of the CH2 decays to CH, which, according 



to IHuebner et al.l ( 1992 ) has a very short lifetime (108 s; rate of 9.3xl0 -3 s -1 ). Formally, a 4- 
generation lHaserl model would best describe the C I emission from CH4. However, the progression 
of long to sh ort pa rent-daughter lifetimes allows us to emulate the distribution with a standard 
2-component lHaserl model, with the effective parent lifetime equal to the sum of all of the lifetimes 
(see Table E]). We mount a similar argument t o estimate limit s on the effective ejection velocity 



of carbon from the chain of CH4 decay. From IHuebner et al.l . the first and last links are 3 and 



0.7kms _1 , respectively, which a dd in quadrature to find 3.1 kms -1 . Excess energies of the CH2 
dissociations are not quoted by Ivan Dishoeck et al.l ( 1990 ). so 3.1 kms -1 is a lower limit to the 
carbon ejection velocity from the CH4 chain. When we created the model described below, we used 
3.7 kms -1 , which assumes a 2 kms -1 ejection velocity of CH from CH2. 

The CH3OH lifetime is 8.7x l0 4 s (l.lxlO -5 s -1 ) and dissociates primarily to H2CO, which 



lives 4600 s (IHuebner et al.lll992l ). H2CO primarily dissociates to CO, so CH3OH cannot provide 
the short-lived component of the C I emission. The ejection velocities of the daughters from the 
CH3OH and H2CO dissociations are 1.7 and lkms" 1 , respectively, so, adding in quadrature with 
the carbon ejection velocity from CO, we find a total v e = 4.5 kms -1 . 

The fact that the final carbon ejection velocities we calculate for our three candidate C I parent 
molecules all lie within 25% of each other is significant. This means that after the details of the 
parent molecule dissociation are sorted out, the carbon from all of these sources will be traveling 
at rou ghly the same speed away from the nucleus of the comet. At large cometocentric distances, 
Haserl models of these individual components will therefore be parallel to each other. This p erhaps 
explains why the radial profiles shown in Figure [3] and Figure [6] follow the characteristic lHaser 
profile outside of ~3xl0 5 km despite the fact that multiple sources of C I are present. 

As a final test of our analyses, we constructed a three-component model consisting of CO 
Fourth Positive and two C I sources, one from CO and one from CH4. The lifetimes and ejection 
velocities are discussed above and the best-fit production rates are shown in Table [5j The grism- 
mode residual is shown in the lower right-hand panel of Figure [5l The model, with its constituent 
parts is shown compared to the direct-mode quadrant 1 profile in Figure El 

A final adjustment to achieve the best residuals to the grism-mode data was made by shifting 
the comet image 2 pixels in the positive dispersion direction. This corresponds to 3", which is 
smaller than the 4.5" GALEX FUV point-spread function and reasonable given the difficulty in 
finding high-precision astrometric solutions in the dispersion direction from stellar spectra at a 
spectroscopic resolving power (A/AA) of 200. 

Having achieved a good residual grism-mode image, we e stablished an upper l imit of Q(S) = 
2xl0 26 s~ 1 , assuming the dom inant parent of S is H2S (e.g., iMcPhate et al.lll999l ) with lifetimes 
from iMeier fc A'Hearnl (jl997h . 
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Fig. 6. — Radial profile from quadrant 1 using 500 counts per bin compared to the three-component 
model discussed in the text. Residuals of this model from the grism-mode data are shown in the 
lower right-hand panel of Figure [5j 
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5. Discussion 

We have created a model of CO and C I emission that describes the GALEX direct and 
grism-mode observations of comet Machholz (Figs. [5H6|). The model shows the surprising result 
that carbon from CH4 dominates carbon emission from CO out to ~2xl0 5 km. At even greater 
cometocentric distances C I emission from CH3OH may be significant. Since the carbon ejection 
velocities from all of these parent molecules is likely to be comparable (Table [5]), the net effect of 
multiple carbon parent molecules on our carbon ionization lifetime measurement is negligible, as 
long as the GALEX FOV samples the region of the coma beyond the effective combined parent 



scale length (indicated by the "knee" in the two-component lHaserl model). Inspection of Figured] 



suggests that such a knee has been reached in the data at ~2xl0 km. 

Carbon emission from CH3OH, which dissociates to carbon through CO, will everywhere have 
a flatter distribution than carbon emission from CO and CH4. If this emission is significant, the 
flat part of the parent distribution may flatten the total C I distribution. If this were the case, the 
carbon ionization lifetime measured by our method would be artificially long by some unknown 
amount. However, we estimate the total C I emission from CH3OH will be no more than that of 
C I from CO (Tabled]), so its contribution will be smaller inside some critical radius. This critical 
radius will likely be larger than the cross-over between the CH4 and CO C I sources. If this cross- 
over occurs outside the GALEX FOV, it will have no effect on our measurements. If the cross-over 
occurs within the GALEX FOV, our estimate of the carbon ionization lifetime could be lowered by 
no more than 12% (the difference between the outflow velocities of carbon from CO and CH3OH). 
This is well within the systematic uncertainty induced by our background subtraction (£ )3.1[ Fig. 
§4.21 and Table [3]) and comparable to the uncertainty in our lifetime measurement caused by our 
choice of r m i n ( §4.ip . 



It is also interesting to note that the CO lifetime may be shorter than 5xl0 5 s (rate > 
2xl0~ 6 s~ 1 ) and th at C I emission fr om the CH4 chain contributes less than we estimate. Using 
the solar maxi mum ITozzi et ajj (119981 ') CO photodissociation rate value and the solar wind ioniza- 
tion rates from lRubin et al.l (|2009l ). the CO lifetime could be as short as 3.3xl0 5 s (3.0xlCT 6 s _1 ). 



When we make this change to our model we find that carbon production from CH4 decreases to 
1.0xl0 28 s _1 but that there is no change in the carbon production rate from CO. The production 
rate of carbon from CO could be lowered toward the Q(C) max value quoted in Table [5] if we in- 
cluded carbon from CH3OH in our model. As we have shown in the previous paragraphs, these 
adjustments do not effect the ultimate goal of this work: the measurement of the carbon ionization 
lifetime. We look forward to continuing the analysis of our grism data to directly measure the 
CO lifetime and better constrain the production rates of the various carbon parent molecules with 
these UV data. Calculation of the excess energies of dissociation of CH2 to solar photodissociation 
would also be useful in verifying our CH4 chain analysis. 
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6. Conclusion 

We have used high-quality FUV observations of comet C/2004 Q2 (Machholz) collected by the 
GALEX satellite and simple comet coma modeling techniques to measure the lifetime of ground 
state atomic carbon against ionization processes in interplanetary space. Scaled to 1 AU, we mea- 
sure a value of 7.1 - 9.6xl0 5 s (ionization rate of 1.0 - 1.4xl0 -6 s _1 ). This measurement compares 
favorably to calculations of the total ionization lifetime from three processes (in order of impor- 
tance): solar wind proton charge exchange, solar photoionization, and solar wind electron impact 
ionization, assuming the comet was in the slow solar wind (Table [3]). 

Comet Machholz was observed at a heliographic lat itude of 30°, which is near the boundary 



between the fast and slow solar wind (jBzowski et al.l l2003). Our measurement was sensitive enough 



to determine that the comet was in the slow solar wind ($5] and Fig. U]). Together with ion tail 
observations, observation of the scale length of C I can therefore provide diagnostics of solar wind 
conditions at heliographic latitudes not often sampled by spacecraft. 

The principal measurement uncertainty in the determination of the carbon ionization lifetime 
was related to the uncertainty in determining the background in the direct-mode comet image. 
The C I scale length was significantly larger than the FOV radius, so the background could not be 
determined a priori from the comet image i t self. A background image recorded a month earlier 



was used, but, as first noted by lSujatha et al.l (120091 ). variable airglow contamination makes simple- 



minded background subtraction of GALEX data recorded at different times highly inaccurate. 
We created a general method for the removal of airglow contamination that accurately considers 
observing geometry. This method could be improved using more physically accurate airglow emitter 
profiles from ionospheric and plasmaspheric models. 

Our method of simultaneously using the GALEX direct and grism-mode observations has made 
it possible to measure the ionization lifetime of carbon and begin to probe the distribution of the 
various carbon parent molecules. With careful cleaning of background sources from the grism data, 
this work will continue. In particular, it will be possible to use our parameter optimizer to directly 
measure the CO lifetime. With the addition of optical depth effects to our ^-factor calculations, 
we will also be able to use the GALEX data to determine if CO comes from a distributed source 
in comet Machholz. 

GALEX has observed two other comets to date in the FUV: 8P/Tuttle and C/2007 N3 (Lulin) 
and six comets in the NUV: C/2004 Q2 (Machholz), 9P/Tempel 1, 73P/Schwassmann-Wachmann 
3 Fragments B and C, 8P/Tuttle and C/2007 N3 (Lulin). The GALEX NUV channel includes 
emission from the bright OH 3080 A band, the CS 2576 A band, dust continuum and CO + features. 
The combined image and grism-mode analysis described here is also applicable to the GALEX 
NUV data. 
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